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CN ■ This paper describes a concept for the high-accuracy absolute calibration 



of the instrumental polarization introduced by the primary mirror of a 
O I large-aperture telescope. This procedure requires a small aperture with 

polarization calibration optics (e.g., mounted on the dome) followed by a lens 
that opens the beam to illuminate the entire surface of the mirror. The Jones 
^. ■ matrix corresponding to this calibration setup (with a diverging incident 

beam) is related to that of the normal observing setup (with a collimated 
^ ' incident beam) by an approximate correction term. Numerical models of 

parabolic on-axis and off-axis mirrors with surface imperfections are used to 
q ! explore its accuracy. © 2008 Optical Society of America 
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c*2 ■ 1. Introduction 

Astronomical polarimetry requires frequent calibration operations to remove the instrumen- 
^ ■ tal polarization introduced by the various optical elements encountered by the light beam 
along its path. This contamination is usually determined by placing calibration optics early 
in the light path, which is used to feed light in a known state of polarization into the instru- 
ment. By measuring the polarization of the light that comes out of the system, it is possible 
to characterize it in terms of its Jones matrix (or, alternatively, the Mueller matrix if one is 
using the Stokes formalism). 

Obviously, we can only characterize and calibrate the optical elements that the beam 
encounters after the polarization calibration optics. Therefore, calibrating a telescope requires 
placing such optics at the telescope entrance, before the first reflection on the primary mirror 
(Ml) occurs. This approach 1 is being successfully employed for the Dunn Solar Telescope 
(at the Sacramento Peak Observatory, managed by the National Solar Observatory) and 
the German VTT on the island of Tenerife (at the Observatorio del Teide of the Instituto 
de Astrofisica de Canarias). In both cases, an array of linear polarizers and retarders are 
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slided in the light path, on top of the entrance window, for calibration operations. Separate 
calibrations are obtained for the telescope and the instrument, so that the former does not 
need to be done as frequently. The Jones matrix of the complete system is then obtained as 
the product of the telescope and instrument matrices. 

Unfortunately, entrance window polarizers are not practical for apertures larger than ~lm 
diameter. In the past this has not been a major concern because: a)solar telescopes have 
apertures that do not exceed 1 m; and b) large astronomical telescopes have not been used 
for polarization measurements. However, this scenario is starting to change. Polarimetry 
is proving to be a very powerful tool to explore a broad range of Astrophysical problems, 2 
resulting in a rapidly increasing interest to develop polarimeters for existing large telescopes. 
Second, the development of the Advanced Technology Solar Telescope (ATST) 3,4 demands 
a reliable method to calibrate a large telescope for high accuracy spectro-polarimetry. 

Solar telescopes have been calibrated in the past by observing magnetic structures and 
making assumptions on the underlying physics. This poses important challenges, however, 
especially when pushing the envolope towards new observational domains. Consider for ex- 
ample the ATST, which is intended to do polarimetry at the ~10 -4 accuracy level. One of the 
common assumptions that is usually made in solar polarimetry is to consider that the con- 
tinuum radiation is unpolarized. However, scattering processes can polarize the continuum 
and generate signals of the order of ~1% in the blue side of the visible spectrum. It has been 
stated 5 that "the direct observation of the polarisation of the continuous radiation still is a 
major outstanding observational challenge". Considerations on the symmetry properties of 
Stokes profiles are not appropriate either. Gradients in the line-of-sight velocity or magnetic 
field introduce spurious asymmetries that can invalidate these assumptions. Furthermore, 
some physical processes operating in the atoms are known to induce Stokes asymmetries 
even in the absence of gradients. Spectral lines forming in the incomplete Paschen-Back 
regime become asymmetric. 6 Moreover, the alignment-to-orientation conversion mechanism 7 
may also cause profile asymmetries. 

In summary, it is very important to have an absolute calibration that does not rely on 
preconceived ideas on the objects under study. This is specially true when the instrumen- 
tation is being used to explore new scientific realms. The present work has been motivated 
by the challenge of calibrating the 4-meter primary mirror of the ATST to meet its very 
stringet polarization requirements. It might also be possible to use the calibration method 
proposed here in other existing large-aperture telescopes. However, the actual design of a 
practical implementation is beyond the scope of this paper. The main point of this work 
is to show that a calibration setup with inclined beam incidence can be used to measure 
the polarization properties of a large-aperture primary mirror. Geometrical effects can be 
calculated and removed from the measured Jones or Mueller matrices resulting in a good 
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approximation to such matrices in normal observing conditions. 

2. The calibration setup 

In this paper we shall consider two different configurations: the normal observing setup (OS) 
and the calibration setup (CS) proposed here. Fig 1 shows a schematic representation of the 
OS and CS for an on-axis Ml mirror. Our ultimate goal is to determine the Jones matrix of 
Ml in the OS (M os ). A direct measurement of M os would require polarization optics of the 
same diameter as the telescope aperture, which is not practical due to technical difficulties. 
The M cs matrix, on the other hand, can be determined by mounting appropriate calibration 
optics (and a mechanical control system) at a height H ca i over Ml (as shown in Fig 1, right). 
A small aperture on the dome is probably an ideal location for it. 

The CS requires some (small) amount of additional optics with respect to the OS. At least 
two lenses are required: one at the entrance to open the beam and another at the detector for 
imaging. These additional elements should be designed with stringent polarization require- 
ments. Aberrations, chromatism and other image imperfections can be tolerated in these 
components, which will only be used for calibration. Even if they introduce some residual 
polarization, this should be easily measurable in the laboratory and removed from the Ml 
Jones matrix. 

The size of the calibration optics affects the accuracy of the calibration. In the limit where it 
fills the entire telescope aperture, the CS and the OS are identical and no correction is needed. 
As the calibration aperture becomes smaller, the incidence angles increase resulting in larger 
corrections (and errors). The simulations presented in this paper consider the pessimistic 
limit in which the calibration optics has a diameter approaching zero (a pinhole aperture). 

3. Basic relations 

For the calculations in this paper it is convenient to use cylindrical coordinates with the 
vertical axis along the propagation direction of the incident light beam. The radial coordinate 
p is the distance from the center of the aperture and is the azimuth angle measured from 
an arbitrary reference. Any given point on the surface of the Ml mirror at coordinates (p, 0) 
is characterized by its complex refraction index N(p, 0). We shall consider here the behavior 
of a monochromatic plane wave. This will allow us to describe the instrumental polarization 
of Ml in terms of its Jones matrix M(p, 0). Appendix 5 gives the Mueller equivalent of the 
most important Jones matrices derived in this work. 

Locally, the behavior of Ml in the vicinity of (p, 0) may be approximated by a reflection on 
a flat mirror of homogeneous refraction index N = n + in. This process adopts a very simple 
form in the reference system of the plane of incidence (the plane formed by the incoming ray 
and the surface normal, see Fig 2). Let us denote the components in the plane of incidence 
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with the subindex p and those perpendicular to it with s. In this frame, the Jones matrix of 
the reflection, J\4 ps , is simply: 8 



M ps 



r p 



r s 

Assuming that the refraction index of air is 1, r p and r s are given by: 



VN 2 - sin 2 9 - N 2 cos 9 
Vn 2 — sin 9 + N 2 cos 9 



V 1 V Dill U IV LVJD U , 

r P = . . ~ , ~ (2) 



cos 9- VN 2 -sin 2 9 
s cos9 + VN 2 -sin 2 9' 
where 9 is the angle of incidence. This matrix can be transformed to the global reference 
frame (p, 0): 

M{p,<t>) = R{-<l>)M P sR{<l>), (4) 
where R(4>) is the usual rotation matrix. Writing down M(p, 0) explicitly: 

/ r p cos 2 + r s sin 2 (r p + r s ) sin0cos0 \ 

= ( / \ • j. x ■ 2 x , 2^ ( 5 ) 

y (r p + r Sj ) sm0cos0 r p sin + r s cos / 

4. Mirror models 

4- A. On- axis mirror 

Consider an axi-symmetric mirror illuminated by a collimated beam (OS). The angle of 
incidence 9 is constant along concentric rings in the mirror (9 = 9{p)). For example, a 
parabolic mirror of focal length F is characterized by the condition: 

tanfl = ^. (6) 

Suppose that the complex refraction index is constant over the surface of the mirror 
(perfect mirror). In this case, r p and r s are only functions of p (because 9 = 9{p)). The 
dependences of M(p, 0) in Eq (5) are easily separable. The Jones matrix of a thin ring of 
radius p is simply: 

M(p) = p M{p, 0)rf0 = pvr f {rS + rp) ° | . (7) 

Jo \ (rs + rp) 

Eq (7) represents the Jones matrix of a non-polarizing system (identical reflectivity and 
retardance for both components of the electric field). This is a well-known property of axi- 
symmetric systems. Notice, however, that the symmetry is broken if one observes away from 
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the center of the field of view. Imperfections in the mirror (irregularities in the refraction 
index caused by coating degradation, dust, etc) may also break the symmetry of the system 
and introduce instrumental polarization. 

Let us now turn to the more general case of an imperfect mirror, defined as one with 
N = N(p,<f)). If this is the case then r p and r s vary across the mirror and Eq (7) is no 
longer valid. We seek to determine a suitable calibration by means of the (measurable) M cs 
matrix. The differences between M os and M cs are due to the different incidence angle of 
the beam (represented in Fig 2). We can expand r p in a power series of a as: 

dv cy 2 d 2 v 

r P (9cs) = r p (9 os - «) = r p (6 os ) ~ a -jjf\o os + Y~d6*\ 0os + " ' ' ^ 
and similarly for r s . Inserting this expansion into Eq (5) we obtain: 



M os (p,0) ~ M cs (p,0) -a 



di iP cos (j) + di, s sm (di, P + c?i, s ) sin cos < 

(di tP + di jS ) sin (j) cos (p di jP sin 2 <fi + di tS cos 2 i 



„2 / A. „„c2 X , A. 



a 2 ( d 2)P cos 2 (j) + ^2,s sin (<i 2 ,p + d 2 ,s) sin</>cos< 

2 I (^2,p + ^2,s) sin cos d 2 , p sin 2 + <i 2iS cos 2 • 



where di tP and di tS have been introduced for notational simplicity: 

d % v 

di, P = ~^Qf\e s j (10) 

(and similarly for In the equations above, a, d itP and dj )S are all functions of (p, </>). 
The angle a can be easily determined from geometrical considerations. However, d^ p and di }S 
are affected by imperfections in the mirror that change over time. Let us separate d^ p into 
two components: a nominal d^ p derived from Eqs (2) and (10) with a theoretical refraction 
index N nom (e.g., from manufacturer specifications), and an unknown 8d^ p due to coating 
degradation, dust accumulation, etc: 

di p d{ p -\- 5di p , (11) 

(and similarly for d iyS ). Inserting this into Eq (9) and integrating over the entire mirror 
surface, we have: 

M os = / / pM os (p, 6)dpd<p = M cs + AM + SM (12) 
Jo Jo 

where p ma x is the radius of the Ml mirror. AM can be calculated numerically as the integral 
of AM{p,9): 



AM(p, 6) = -a 



di p cos 2 (f> + dx >s sin 2 (f> + d liS ) sin cos < 

(di tP + di >s ) sin (j) cos (p d\ jP sin 2 <fi + cos 2 i 
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+ 



a 2 I d,2, p cos 2 <fi + d-2,s sin 2 <fi 
2 I (d 2 , p + d 2 , s ) sin cos 



(d 2jP + d 2)S ) sin <f) cos 
o?2,p sin 2 cj) + d,2 7 s cos 2 



) 



+ 



(13) 



Mirror imperfections are accounted for by the (measured) M cs , whereas non-collimated 
incidence is accounted for by the (calculated) AM. 5M is an unknown second-order term 
that couples mirror imperfections and non-collimated incidence. This term is small (as shown 
below) and may be neglected for our purposes here. 

The number of terms to retain in the Taylor expansion of Eq (13) depends on the particular 
telescope configuration and the accuracy required. Typical examples are presented below in 
which AM can be neglected entirely (on-axis mirror) or needs to be calculated up to second 
order (off-axis mirror, see §4.B). 

In the reminder of this section I present the results of numerical simulations that provide 
some insight into the various terms that are involved in the calibration procedure. The 
parameters of the simulation are listed in Table 1. They represent a 4-m on-axis telescope 
with a silver coating on the Ml mirror. The coating has been degraded so that the complex 
refraction index fluctuates over the mirror surface as N(p, <fi) = N av [l + N^cos(4>/2)} (N av is 
the average refraction index and N* is the amplitude of the fluctuation). This choice has been 
made to represent a pessimistic scenario that induces a considerable amount of instrumental 
polarization. The discretization of the simulation considers 100 points in p and 200 in <p. 

The Jones matrices M os and M cs , obtained by applying Eq (5) to each area element of 
the mirror, are: 



The "irregularities" introduced in the refraction index of the mirror break the symmetry 
and give rise to polarizing effects in M os , with off-diagonal terms of approximately 6% . 
Fortunately, the calibration setup matrix M cs is an excellent approximation to M os , with 
a maximum difference of ~10 -5 . It is then possible to calibrate the telescope almost down 
to the 10~ 5 level without even having to correct for the inclined incidence (i.e., neglecting 
AM in Eq [12]) 

4-B. Off- axis mirror 

An off-axis mirror can be represented by a larger "equivalent" on-axis mirror with a variable 
refraction index, as depicted in Fig 3. The equations derived in §4. A above are still valid in 




(14) 




(15) 
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the (p, 0) reference frame. One simply needs to set the reflectivity to zero outside the shaded 
area of the figure (r > r max ). This can be accomplished, e.g. by setting the refraction index 
to 1. 

This type of mirrors is slightly more complicated to calibrate. Even a perfectly coated 
mirror will produce instrumental polarization. A simulation with the parameters listed in 
Table 2 yields the following matrices: 

M os = -0.97exp(0.59i) 

M cs = M os + L24 x io-2 exp ( 2 .22z) [ L0 ° ) . (17) 

y 0.00 0.95exp(3.11z) J 

M cs is now significantly different from M os and we need to calculate the correction term 
AM (Eqs [12] and [13]). Let us consider for the moment that the refraction index is perfectly 
known (5di iP = 5di tS = in Eq [11]). Calculating AM with Eq (13) up to second order we 
have that: 

M cs + AM = M os + 2.22 x 10" 4 exp(2.65z) | 1,0 ° |. (18) 

V ; ^ 0.00 0.43exp(0.44z) J v ; 

Let us now turn to the more general case of an off-axis mirror with surface irregularities 
for which we have only an imperfect knowledge of the average refracion index. Again, we use 
a refracion index with an angular dependence N(p, <ft) = N av [l + N*cos(<f)/2)]. Furthermore, 
we do not know exactly the average refraction index of the mirror N av , but only an approx- 
imation N nom . This approximate value will be used in the calculation of di jP , di )S and AM 
(Eqs [11] and [13]). I carried out several experiments with different values of N nom and 
to determine the sensitivity of the calibration to these parameters. Some of the calculations 
include only the first-order dependence of AM on a (first term in the right-hand side of 
Eq [13]). These are denoted by AM* (as opposed to AM, which considers the second-order 
term). Table 3 lists the results of these experiments. The first two rows give the highest 
polarizing term in M os for each simulation. The third to fifth rows show how the calibration 
error decreases with successive levels of approximation. By comparing the second and fifth 
rows, we can see that the calibration is able to reduce the instrumental polarization by an 
amount between one and two orders of magnitude. 

5. Conclusions 

This paper introduces a new concept to calibrate telescopes for astronomical polarimetry. 
The proposed method is particularly useful for modern large-aperture telescopes, for which 



0.00 
0.00 



0.00 

-3.51 x 10" 2 exp(1.61z) 



(16) 
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it is probably the only practical procedure (at least for purely instrumental calibration). 
An accurate absolute calibration will be crucial for the new weak-signal science that the 
ATST will open. Existing night-time telescopes may also take advantage of this calibration 
procedure. 

The Jones (and Mueller) matrix of an on-axis mirror is almost unaffected by the non- 
collimated incidence of the beam in the CS. For the particular configuration considered in 
§4. A, the Jones matrix obtained from the calibration is good to almost 10~ 5 with no need 
for any additional correction (see Eq [15]). This is in spite of the relatively large mirror 
imperfections in the simulation, which generate polarizing terms of the order of 6%. 

An off-axis mirror suffers more instrumental polarization due to the asymmetric configura- 
tion. A mirror that produces instrumental polarization of a few percent can be calibrated to 
reach the 10~ 4 level (see Table 3). In addition to measuring the Jones matrix of the calibra- 
tion setup (M cs ), it is also necessary to calculate the correction term AM. This calculation 
is straightforward, though, and AM does not need to be recalculated unless the mirror is 
recoated or it degrades to a point where its average refraction index changes significantly. 
Note that the 10~ 4 calibration accuracy includes some uncertainty on the average refraction 
index of the mirror. 

Appendix A: Mueller formalism 

I have used in this paper the Jones matrix formalism, which deals directly with the com- 
ponents of the electric field of the light wave. Sometimes, however, the Mueller formalism 
is more adequate, especially when dealing with partially polarized or non-monochromatic 
light. Many researchers are more familiar with the Mueller matrices and the Stokes param- 
eters. For these reasons it is probably useful to provide the Mueller equivalent 9 of the Jones 
matrices derived in this work. 
Eq (14) becomes: 



0.89 
-3.42 x 10~ 5 
1.65 x 1(T 2 



-3.42 x 1(T 5 

0.88 
5.49 x 10~ 7 



1.65 x 10^ 2 
5.49 x 10- 7 



2.19 x 10^ 6 \ 

0.11 
1.45 x 10" 5 



AT 



OS 



(Al) 



V 2.19 x 10 
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0.11 



0.89 
1.45 x 10~ 5 



0.88 



/ 



Eq (15) becomes: 
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M 



cs 



M 



OS 



-2.23x10" 



0.30 

1.00 
-6.77 x 10" 2 
-6.42 x 10~ 2 



1.00 

0.30 
-1.80 x 10~ 2 
-5.45 x 10~ 2 



-6.77 x 10~ 2 
-1.80 x 10" 2 
0.30 
-0.45 



6.42 x 10~ 2 \ 



5.45 x 10 
0.45 
0.30 



-2 



Eq (16) becomes: 



M 



OS 



( 0.94 
-1.93 x 10~ 3 
2.29 x 10~ 7 
0.00 



-1.93 x 10" 
0.94 
0.00 
0.00 



2.29 x 10~ 7 
0.00 
0.94 
-3.29 x 10~ 2 



0.00 
0.00 
3.29 x 10~ 2 
0.94 



\ 



Eq (17) becomes: 



M 



cs 



M os + 2.34 x 10" 



x 



j 6.92 x 10~ 3 5.85 x 10~ 2 0.00 

5.85 x 10~ 2 6.92 x 10~ 3 0.00 

0.00 0.00 2.96 x 10~ 2 

V 0.00 0.00 1.00 



(A2) 



(A3) 



0.00 
0.00 
-1.00 
2.96 x 10~ 2 



Eq (18) becomes: 



(A4) 



M CS +AM 



M O5 +1.82xl0" 4 x 



-0.97 
-0.13 
-1.78 x 10~ 5 
0.00 



Acknowledgments 



-0.13 
-0.97 
0.00 
-1.13 x 10~ 5 



-1.78 x 10~ 5 
0.00 
-1.00 
-0.72 



0.00 \ 
1.13 x 10" 5 
0.72 
-1.00 J 

(A5) 
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Table 1. Simulation of a 4-m on-axis mirror 



Parameter Value 

Pmax 2 771 

F 8 m 

H ca i 10 m 

N av 0.2 + 3.4z 

N* 0.50 
max(|M os - M cs \) 1.56 x 10~ J 
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Table 2. Simulation of a 4-m off-axis mirror 



Parameter Value 

Tmax 2 771 

Pmax 6 m 

F 8 m 

H ca i 10 m 

N av 0.2 + 3.4i 

N* 0.0 
max(|M° 5 -M cs |) 1.24 x 10" 2 



12 



Table 3. Simulation of a 4-m off- axis mirror with surface irregularities 



Nf = 0.10 Nf = 0.10 N* = 0.25 = 0.25 

J^nom _ pjav jynom _ \ Q\N av ]\[ nom = l.OliV™ N n ° m = 1.10iV a 

\ M °H ~ \ M ?.i\ 2.02 x 10~ 3 2.02 x 10~ 3 2.25 x 10~ 3 2.25 x 10" 3 

|Mgf | 1.50 x 10" 2 1.50 x 10" 2 3.86 x 10~ 2 3.86 x 10~ 2 

max{\M os - M cs \) 1.25 x 10~ 2 1.25 x 10" 2 1.26 x 10~ 2 1.26 x 10~ 2 
max[\M os - 

(M cs + AM*)|] 4.04 x 10~ 3 3.89 x 10~ 3 3.71 x 10" 3 2.59 x 10~ 3 
max[\M os - 

(M CS + AM)\] 2.84 xlO" 4 4.18 x 10" 4 7.11 x 10~ 4 1.80 x 10~ 3 
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Fig. 1. Left: Normal observing configuration (OS), with a collimated incident 
beam. Right: Calibration configuration (OS), with an inclined (diverging) in- 
cident beam, navarrofl.eps. 



Fig. 2. Incidence angles 9 and a for the OS and CS, respectively. Solid (dashed) 
represent rays in the OS (CS). The dotted line represents the normal to the 
mirror surface, navarrof2.eps. 



Fig. 3. Schematic representation of the equivalent on-axis mirror. Left: Lateral 
view. The thick line represents the actual off-axis mirror of radius r max . The 
thin line represents the equivalent on-axis mirror. Right: Top-down view. The 
shaded area is the actual off-axis mirror, navarrof3.eps. 
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